Mesoscopic effects in adiabatic spin pumping 
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We show that temporal shape modulations (pumping) of a quantum dot in the presence of spin- 
orbital coupling lead to a finite dc spin current. Depending on the strength of the spin-orbit coupling, 
the spin current is polarized perpendicular to the plane of the two-dimensional electron gas, or has 
an arbitrary direction subject to mesoscopic fluctuations. We analyze the statistics of the spin and 
charge currents in the adiabatic limit for the full cross-over from weak to strong spin-orbit coupling. 
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There is a growing interest in the physics of spin trans- 
port through low-dimensional quantum structures [1] 
with the aim of controlling and manipulating spin in mi- 
croelectronic devices. A fascinating tool for spin manip- 
ulation is a "spin battery" or "spin pump" , a device that 
generates a spin current without an accompanying charge 
current [2] . Following a proposal of Mucciolo et al. [3] , a 
"spin battery" was realized recently [4] using a quantum 
dot in a two-dimensional electron gas (2DEG). In this 
geometry, the direction of spin polarization is set by an 
external magnetic field parallel to the plane of the 2DEG. 
Current is generated by periodic variation of gate volt- 
ages, an eventual charge current being suppressed by fine 
tuning the dot shape. Other proposals for spin pumps 
have used the idea of locally breaking spin-rotation sym- 
metry to pump a spin-polarized current by a magnetic 
field [2], magnetic impurities [5, 6, 7], and by periodic 
modulation of the spin-orbit coupling in the 2DEG [8] . 

Spin current is different from charge current because of 
the vector nature of spin. In fact, it is the "vector" nature 
of spin — different directions of spin corresponding to dif- 
ferent quantum-mechanical superpositions of "spin up" 
and "spin down" — that makes it a promising candidate 
for the practical realization of a quantum computer [9]. 
Control of magnitude as well as direction of spin current 
is of paramount importance if the full benefits of a spin- 
tronic circuit are to be reaped, or if the electron spin is 
used as the building block in a scheme that makes es- 
sential use of quantum coherence. The "spin battery" of 
Rcfs. 3, 4 satisfies this requirement in part, since mag- 
nitude and sign of the current can be controlled, but its 
direction is determined by the external magnetic field, 
allowing for spin-polarization in 2DEG plane only. 

In this letter, we investigate a quantum-dot based 
spin pump for which Bychkov-Rashba [10] and Dres- 
selhaus [11] contributions to spin-orbit coupling are the 
sources for spin-rotation symmetry breaking. We show 
that, depending on the strength of the spin-orbit cou- 
pling, a quantum dot with a variable shape can pump 
either a spin current perpendicular to the 2DEG plane, 
with the sign and magnitude of the current subject to 
control via mesoscopic fluctuations, or a spin current with 
an arbitrary direction. 



The possibility to generate a spin current perpendic- 
ular to the plane of the 2DEG allows for an interesting 
realization of a "spin Hall effect" [12] if the spin current 
is injected into a 2DEG which shows the anomalous Hall 
effect [13] . The 2DEG then has skew spin-orbit scattering 
with spin up and down (measured perpendicular to the 
2DEG plane) being scattered preferentially in opposite 
directions, perpendicular to the flow of current. While a 
current of unpolarized electrons leads to a spin imbalance 
perpendicular to the current flow [12] and a spin current 
polarized in the plane of the 2DEG is not affected, a spin 
current polarized perpendicular to the 2DEG should lead 
to an anomalous Hall voltage across the sample. 

The system under consideration consists of a ballis- 
tic quantum dot connected to two electron reservoirs 
through ballistic point contacts with TVi and N2 channels 
each, see Fig. 1, inset. Gate voltages xi and X2 of two 
shape-distorting gates allow for a time-dependent varia- 
tion of the dot shape. The same geometry was used in 
the "spin battery" of Ref. 4, with a quantum dot without 
notable spin-orbit coupling. Electron motion in the dot 
is characterized by the transit time rtr = L/vp, where 
L is the dot size and vp the Fermi velocity, and the 
time Tcsc = h/NA for escape to the reservoirs, where 
N = Ni + N2 and A is the mean spacing between single- 
electron levels in the quantum dot. The escape time, 
which is determined by the point contacts, is typically 
much larger than Ttr. In the absence of spin-orbit cou- 
pling, periodic variation of the gate voltages xi and X2 
leads to a dc charge current through the dot [14, 15, 16]. 
Below, we show that a spin current is generated in the 
presence of spin-orbit scattering. 

The Dresselhaus and Bychkov-Rashba contributions to 
the spin-orbit coupling arise from the spin splitting of the 
conduction band in bulk GaAs [17] and the asymmetry of 
the potential well forming the 2DEG, respectively. The 
corresponding term in the Hamiltonian has the form 



■f^SO = Q [PxO-x - PyCTy] + a [p X a]. 
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where px and py is the in-plane electron momentum and 
ax, (Ty, and are the Pauli matrices. The coefficients g 
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and a define a length scale 



A 
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where m is the electron mass. For quantum wires, only 
the spin projection perpendicular to the wire and in the 
plane of the 2DEG is conserved. As a result, pumping 
through such a system yields a spin current with the po- 
larization direction in the plane of the 2DEG [8]. For 
quantum dots of size L X, the role of spin-orbit scat- 
tering is qualitatively different from one-dimensional or 
bulk systems. The difference originates from the fact 
that, although spin-orbit scattering has a small effect in 
the transit time rtr if i <C A, it still may have a large 
effect within the time scale Tosc [18]. As shown in Ref. 
18, a unitary transformation casts H into the form 



(2) 



where 



^ , a|| = —{ar-a + g[r x a\^} (3) 

are spin-dependent vector potentials representing the ef- 
fects of spin-orbit scattering. They are characterized by 
scattering times tj_ and r|| , 

Ti_ = /^-Vtr(2A/L)^ r|| = (K')"'^i(2A/L)' » r^, (4) 

where n and k' are geometry-dependent coefficients of 
order unity [19]. Whether or not the spin-orbit term 
ay is important depends on the relative size of r|| and 
Tosc- We now discuss the cases ry ^ Tosc and T|| < Tcsc 
separately. 

For r|| 3> Tesc the only relevant spin-orbit term in the 
Hamiltonian is the spin-dependent vector potential aj_. 
This term has the same form as the vector potential aris- 
ing from a magnetic field of size Bso = /i/2eA^ with op- 
posite directions for spin up (t) and spin down (|), mea- 
sured perpendicular to the plane of the 2DEG. Hence, the 
perpendicular spin projections are conserved, and sepa- 
rate dc currents and /| will be pumped through the 
quantum dot for spin up and spin down. Generically, 
one has a spin current /^ ^ — Iij ~ 7^ because electron 
trajectories that enclose the same geometrical area (see 
inset of Fig. 1) gather an "Aharonov-Bohm" flux with op- 
posite signs for electrons with opposite spin projections, 
thereby giving the scattering phase shift a spin depen- 
dence. In a chaotic cavity, both the sign and magnitude 
of the currents I^^i is essentially random, dependent on 
the detailed shape and electron density of the quantum 
dot. Following the protocol of Ref. 3, a spin battery with 
spin current perpendicular to the 2DEG plane results if 
the shape of the dot is fine tuned so that the charge cur- 
rent /c = /| + /| = 0. In order to compare typical (rms) 
sizes of spin and charge currents, we note that the spin 



current is formally equal to the magnetic-field antisym- 
metric component of the charge current for a quantum 
pump with spinless electrons in a perpendicular magnetic 
field of size Bgo- The latter quantity was calculated in 
Ref. [20] using random matrix theory, which is valid if 
the dot has an irregular shape and Tosc ^ Ttr, 
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We now turn to the general case Ty < Tosc- We con- 
sider a sinusoidal time dependence of the gate voltages, 
xi(t) = Sxi cos{ujt) and X2(t) = 8x2 cos(a;t + (p) and con- 
sider the contribution to the charge and spin current that 
is bilinear in the gate voltage amplitudes 5xi and 5x2- 
Starting point of our calculation is the relation between 
the charge and spin currents and the 2N x 2N scattering 
matrix S of the quantum dot [14], 
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, (6) 



where A is an A'^ x A'^ diagonal matrix with elements Kjj = 
N2/N for j < Ni, and Ajj = -NJN for Ni < j < N 
and 1 is the 2x2 unit matrix in spin space. 

In order to calculate the average and root-mean-square 
(rms) of the charge and spin current for an ensemble of 
quantum dots, we need to know the average of a product 
of up to four scattering matrix elements taken at differ- 
ent values of the parameters xi and X2- We calculate 
this correlator using random matrix theory (RMT) [21]. 
In RMT, the spin-orbit part of the Hamiltonian (2) is 
replaced by a 2M x 2M random hermitian matrix [19] 



A3 (g) CT3 Ai ® ai + A2 (E) (72 



, (7) 



whereas the dependence on the shape-distorting gate 
voltages xi and X2 is represented through the random 
hermitian matrix 



A' 



shape 



XjXj (g) 1. 



(8) 



Here, Xj, j = 1,2 are real symmetric random M x M ma- 
trices and Aj, j = 1, 2, 3 are real antisymmetric random 
M X M matrices with 



tr = ArSi. 



(9) 



At the end of the calculation, the limit M — > 00 needs 
to be taken. Performing the random matrix average us- 
ing standard methods [19, 22], we find that the relevant 
average of a product of four scattering matrices is 



{S{l)ki<E>S\2'U<S>S{l%, 
= W^^'>{12')W^^\l'2)6in 



.<S)S\2)i^) 
-T^(2)(i2'l'2), 



(10) 
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where "1" , "2" , "1"' , and "2"' are shorthand notations 
for values of the gate voltages xi and X2 and the Fermi 
energy e and we have used tensor notation for the spin 
degrees of freedom. The first contribution on the r.h.s. 
of Eq. (10) is the product of pair averages [19], 



M^(i)(12') 



1 



MI2 -tri?(l) (g)i?t(2') 



12 



with a similar definition of VF(^'(1'2). Here I2 = 1 1, 
the auxiliary matrix R is defined as 

R{e,xi,X2) = cxp [2ni{e - i?so - ^^shapo(a;i, 2:2))] , 

(11) 

and the trace "tr" is not taken over the spin degrees of 
freedom. In taking the inverse in Eq. (11), the rule for 
multiplication of the tensor products is (aj ak> ) (cTfe (8) 
"'i') ^ i^j'^k) ® (cj'CTfc/). The second contribution in Eq. 
(10) involves the random matrix equivalent of the Hikami 
box from diagrammatic perturbation theory and can be 
calculated using the method of Ref. 22. Using the tensor 
notation with the same product rules as before, the result 
can be written as 



W(2)(i2'i'2) = 



l^(i)(12')P^(i)(l'2)P(12'l'2) 
X M^(i'(12)W^(i)(l'2'), 



(12) 



where 



I?(12'1'2) = MI4 - tr i?(l) R\2') ® R{1') (g) R''{2), 

The traces are calculated using Eq. (9) after expanding 
R to second order in Hgo and -ffshapc and taking Af — *■ 00. 

Substituting these results into Eq. (6) one then finds 
the average and variances of the spin and charge currents. 
All ensemble averages are zero, as well as the covariances 
for different components of the spin current or of spin and 
charge currents. The full results for the average square 
spin and charge currents can be written by introducing 
Jo = {uj/2Tr)4:Sxi6x2 sm(f>y^¥hN^/I^, c± = Tcsc/t±, and 



{ID 
(IL) 
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(l + 2c||)3 (1 + C||+Ci)3 
1 2 



° ' (1 + 2C||)3 (1 + C|| +Ci)3^ 



These general results confirm that spin current is po- 
larized perpendicular to the plane of the 2DEG if r|| ^ 
Tosc, while its direction is arbitrary if r|| <C Tesc- In prac- 
tice, the ratio T\\/Tesc can be tuned rather straightfor- 
wardly by changing the conductances of the point con- 
tacts connecting the quantum dot to the outside world or 
the dot size. Recent experiments by Zumbiihl et al. show 



that both limits can be obtained experimentally [23] : De- 
pending on the size of the quantum dots, Zumbiihl et al. 
find c^N ranging from 0.94 to 20 and c\\N ranging from 
0.25 to 20. For iV = 2 and at a temperature T < t~1, 
we obtain rms /s,z/rms Ic = 0.72. At finite tempera- 
ture, both the effect of thermal smearing and a finite 
dephasing time need to be taken into account. De- 
phasing is accounted for by the substitution 1/tosc — >■ 
l/''csc + l/''"0j whereas thermal smearing requires inte- 
gration over e. For large temperatures T ^ t~^1, and 
in the absence of any Zeeman coupling, we can borrow 
the results of Ref. [20] and find the polarization ratio rms 
/.^^/rms/c = [(ci + 2cA_)l{c\ + 2c^ + 2)]i/2 = 0.61. 

Before concluding, we discuss the dependence of the 
pumped current on an applied magnetic field B. This 
question is relevant if a spin pump is used in conjunction 
with spin manipulation by means of in-plane or perpen- 
dicular fields. Within RMT, the effect of a magnetic field 
is modeled by a third random matrix [18, 19] 



1 




— B-a, (14) 
2rz 



where 1/tz = fisgB, B is the magnitude of the applied 
magnetic field, B its direction, /is is the Bohr magne- 
ton, g is the g factor of GaAs, th ~ K~^Ttr(2/eBzL^)^, 
th,± = {AX/ ^BgLB)^ / K"Tti-, K and k" being geometry- 
dependent coefficients of order unity, ^ is an M x M real 
antisymmetric matrix with trace tr^^""" = M^, and X 
is a real symmetric matric with trace ti X'^ ~ M^. The 
time Th describes the orbital effect of the magnetic field 
component perpendicular to the 2DEG; it is the time 
needed for picking up a quantum of magnetic flux. The 
time th,i_ describes spin-flip processes arising from the 
interplay of the component of the magnetic fleld parallel 
to the 2DEG and the spin-orbit scattering [24]. 

Adding the random matrix Hb to the exponent in 
Eq. (11) and repeating the previous calculations, we flnd 
(i) inclusion of th has no effect on the current statis- 
tics and (ii) inclusion of th,i_ amounts to the subsitution 
1/t± 1/t± + 1/th,±- In general the effect of Zeeman 
coupling to the parallel fleld (time scale tz) is to de- 
crease correlations between J-|- and so that both rms 
Is^z and rms Ic are reduced. On the other hand, the vari- 
ance of the spin current polarized along the direction of 
the in-plane magnetic fleld is enhanced. Both these ef- 
fects arc shown in Fig. 1. The full results for the current 
variance in the presence of Zeeman coupling are rather 
lengthy and will be reported elsewhere. Here, we con- 
flne ourselves to the limiting cases of Ref. [18] which are 
distinguished by a parameter £ = 1,2 characterizing the 
mixing of states with different spins for strong Zeeman 
splitting. Choosing the magnetic fleld along the x direc- 
tion, these are: (i) C|| = Tosc/t|| <C 1, and b = Tcsc/tz 
large or small in comparison to cj_ ~ Tcsc/t±, for which 
S = 1; and (ii) c± ':$> I, b large or small in compar- 
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FIG. 1: Current variance (in units of /o/2) as a function 
of Zeeman coupling strength (6 = Tcsc/tz) for a fixed value 
of SO coupling strengths c± — 0.47, and cy = 0.12 taken 
from Ref. [23]. The difference between {Ib,x) and {Is,z) for 
large b vanishes for 2> 1. Inset shows a quantum dot 
with out of phase shape modulation, and a typical electron 
trajectory which gathers a Aharonov-Bohm flux depending 
on its perpendicular spin projection. 



ison to C||; for which S = 2. For the case (i), when 
b'^ /c± ^ 1 <C cj^; we obtain to the lowest order: 



and when cj^ <C 1 <C & we find: 



(^c ) = (IL) = ^0 (2 - 3c^ + 6ci) 



il 



s.z/ 



3c I — 6c I 



(15) 



(16) 



which is "dual" to expression (15) under the exchange 
Is.x Is,z and cj^ ^ AlP'/c^. For the case (ii), the 
lowest order expansion in b"^ /c±^ <C 1 <C cj_, C|| 



2 



II (1 ± 3/8 cl^ ^ l/2(6Vcx)c 



{I'D = {ID = il (1 - 1/8 c 



-3^ 



(17) 



is dual under the transformation 46^/c^ C|| to the 
variance in the other limit C|| <C 1 <C c^,b'^/c^. The 
non- vanishing current variances (/^) can be written in a 
unified way as (/^) = 4(s//3I])/q, where /3 = 2 defines 
the time reversal symmetry of orbital motion, and s = 1 
is the Kramers degeneracy parameter. 

In conclusion, we have shown that a quantum pump 
consisting of a quantum dot with spin-orbit coupling 
allows for the generation of both a dc charge current 
and a dc spin current. Depending on the strength of 
the spin-orbit scattering or on the conductances of the 



point contacts between the dot and the reservoirs, the 
pumped spin current is perpendicular to the 2DEG or 
has an arbitrary direction. Measurement of the spin cur- 
rent is possible via the spin Hall effect or by connec- 
tion of the spin pump to a (semiconducting) ferromagnet 
with known magnetization direction. Furthermore, the 
method of adiabatic spin pumping, as opposed to other 
mechanisms such as rectification [4, 25], allows the direc- 
tion of spin polarization to be continuously changed. 
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